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1 Introduction 



We consider exactly solvable models of statistical mechanics in one space and one time dimension. 
The Quantum Inverse Scattering Method and Algebraic Bethe Ansatz are effective methods for a 
description of the spectrum of these models. Our aim is the evaluation of correlation functions of 
exactly solvable models. Our approach is based on the determinant representation for correlation 
functions. It consists of a few steps: first the correlation function is represented as a determinant of 
a Fredholm integral operator, second — the Fredholm integral operator is described by a classical 
completely integrable equation, third — the classical completely integrable equation is solved by 
means of the Riemann-Hilbert problem. This permits us to evaluate the long distance and large 
time asymptotics of the correlation function. The method is described in Q, ||. 

The most interesting correlation functions are time dependent correlation functions. The deter- 
minant representation for time dependent correlation functions was known only for the impenetrable 
Bose gas (the spectrum of the Hamiltonian of this model is equivalent to free fermions). In this paper 
we have found the determinant representation for the time dependent correlation function of local 
fields of the penetrable Bose gas. The main idea for the construction of the determinant representa- 
tion is the following. We introduce auxiliary Bose fields (acting in the canonical Fock space) in order 
to remove the two body scattering matrix and to reduce the model to the free fermionic case. We 
want to emphasize that all dual fields, which we introduce commute (belong to the same Abelian 
sub-algebra). Therefore we do not have any ordering problem. This will also permit us to perform 
nonperturbative calculations, which are necessary for the derivation of the integrable equation for 
the correlation function. 

First we shall discuss our model. 

Quantum nonlinear Schrodinger equation (equivalent to Bose gas with delta-function interaction) 
can be described by the canonical Bose fields tp(x) and tp^x) with the commutation relations: 

ty(x),^(y)] = 6(x-y), [^{x),^y)\ = [^(x),^{y)) = 0, (1.1) 

acting in the Fock space. Fock vacuum |0) and dual vector (0| are important. They are defined by 
the relations 

t/j(x)\0)=0, (0|V> f (x) = 0, (0 | 0) = 1. (1.2) 
The Hamiltonian of the model is 

H= dx (j) x ijj^ (x)d x ip(x) + ap^ (x)ip^ (x)ip(x)ip(x) — hip^ (x)ijj(xfj , (1-3) 
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Here c is the coupling constant and h > is the chemical potential. We shall consider the repulsive 
case < c < oo. 

The spectrum of the model was first described by E. H. Lieb and W. Liniger Q, |J|. The Lax 
representation for the corresponding classical equation of motion 

i— 4 = [^H\ = --^ + 2ctfW>-W, (1.4) 
at dx 

was found by V. E. Zakharov and A. B. Shabat ||. The Quantum Inverse Scattering Method for 
the model was formulated by L. D. Faddeev and E. K. Sklyanin 

In this paper we shall follow the notations of 0]. First the model is considered in a finite periodic 
box of length L. Later the thermodynamic limit is considered when the length of the box L and the 
number of particles in the ground state go to infinity, with the ratio N/L held fixed. 

The Quantum nonlinear Schrodinger equation is equivalent to the Bose gas with delta-function 
interaction. In the sector with N particles the Hamiltonian of Bose gas is given by 

N Q2 

n N = -Y.Jf2 +2 ° £ 8(z k -Zj)-Nh. (1.5) 

j=l OZ j l<j<k<N 

Now a few words about the organization of the paper. 

In Section 2 we shall review the Algebraic Bethe Ansatz and collect all the known facts necessary 
for further calculations. In Section 3 we shall calculate the form factor of the local field in finite 
volume. In Section 4 we shall present the idea of summation with respect to all intermediate states. 
In Section 5 we introduce an auxiliary Bosonic Fock space and auxiliary Bose fields. This helps us 
to represent the correlation function as a determinant in the finite volume. In Section 6 we consider 
the thermodynamic limit of the determinant representation for correlation function. Length of the 
periodic box L and number of particles in the ground state go to infinity but their ratio remains 
fixed. This leads us to the main result of the paper (see formulae ( |6.24 )— (|6.27D ). The correlation 
function of local fields in the infinite volume is represented as a determinant of a Fredholm integral 
operator. For evaluation of the thermodynamic limit it is necessary to sum up singular expressions. 
Appendix A is devoted to these summations. In Appendix B we present realization of quantum 
dual fields as linear combinations of the canonical Bose fields. Appendix C shows how to reduce the 
number of dual fields. Appendix D contains determinant representation for temperature correlation 
function. 

In forthcoming publications we shall use the determinant representation for the derivation of 
completely integrable equation for correlation functions. Later we shall solve this equation by means 
of the Riemann-Hilbert problem and evaluate the long-distance asymptotic. 
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2 Algebraic Bethe Ansatz 



Let us review some main features of the Algebraic Bethe Ansatz, which we shall use later. We 
consider the quantum nonlinear Schrodinger model. The starting point and central object of the 
Quantum Inverse Scattering Method is the R-matrix, which is a solution of the Yang-Baxter equa- 
tion. For the case of the quantum nonlinear Schrodinger equation, it is of the form : 



( /OM) 






where 



ic 



X — \i 





1 



/(A,m) 





1 





\ 





(2.1) 



A — fi + ic 
X — \i 



Later we shall also use functions 



t(x, n) 



ic 



3(A,ju) 



(2.2) 



(2.3) 



ic (A — ( A — n + ic) h(X, /j,) 

Another important object is the monodromy matrix 

A(X) B(X) 
C(A) D(X) 

The operators A, B, C, D are acting in the Fock space where the operator ij)(x) was defined. Their 
commutation relations are given by 



T(A) 



(2.4) 



R(X, m) (T(A) T(/x)) = ® T(A)) fl(A, /i). 

These relations are written out explicitly in Section VII. 1 of ||. 
The hermiticity properties of T(A) are 

a x T*(X)a x = T(X), 



(2.5) 



(2.6) 



so that fit (A) = C(A). 

The Hamiltonian of the model can be expressed in terms of A(X) + .D(A) by means of trace 
identities (Section VI. 3 of Q). The vacuum is eigenvector of the diagonal elements of T(X) 



A(A)|0)=a(A)|0) 

<0|A(A) = a(A)(0| 
f iLA 

a(A) = exp 



D(A)|0) = d(A)|0) 
(0|L»(A) = d(A)(0|; 
d(X) = exp j^-J 



(2.7) 
(2.8) 
(2.9) 
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Later we shall also use the function 

KA) = $!=e- (2.10) 

The operator C(A) annihilates the vacuum vector and the operator -B(A) annihilates the dual vacuum: 

C(A)|0)=0, (015(A) =0. (2.11) 

The Hamiltonian of the model commutes with A(X) + D(X) and they can be diagonalized simul- 
taneously. The eigenvectors of the Hamiltonian are 

N N 

Y[B(nj)\0), and (0|nC(^), (2.12) 
3=1 3=1 

if iXj satisfy Bethe Equations 

^frfe^ = 1) or ^fi^Bk = irl) N-i (2 . 13) 



It is convenient to rewrite ( 2.13| ) in logarithmic form. For the ground state 



W+ . = ^ + £«>.(£±^)-*r(,-^) (2.14) 



It is proven in Section 1.2 of || that solutions /Xj of equation (2.14) are real. 

The distribution of \ij in the ground state in thermodynamic limit can be described by linear 
integral equation. The thermodynamic limit is defined in the following way: N — > oo, L — > oo and 
N/L = D is fixed. In this limit fj,j condense (fMj+i — [J-j = 0(1/ L)) and fill the symmetric interval 
[ — where q is the value of spectral parameter on the Fermi surface. 

In the thermodynamic limit the function of local density p(p) can be defined in the following 
way 

p( Mj .)=lim — 1 (2.15) 
L \H+i ~ H> 

The lim in the r.h.s. denotes the thermodynamic limit. This function satisfies the Lieb-Liniger 
integral equation 

P(P) " ^ / K(u, ii)p(v)dv = — . (2.16) 

-9 

Here 

2c 

C 2 + (pi — V 



K{v,n)= „ 2 , - „, 2 , (2.17) 
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and 



N f 
D = — = J dnp{n). 



(2.18) 



In such a way we have described the ground state. 

Now we can define the correlation function of the local fields 



(V'(0,0)V' t (2;,t)) = lim 



N N 

(0\l[C(^(0,0)^(x,t)l[B( H )\0) 
i=i i=± 

N N 

<oi n c <pj) n 

3=1 3=1 



(2.19) 



Here 



tf{x,t) = e iHt ^(x,0)t 



iHt 



(2.20) 



We shall use the notation fj,j for the ground state only. The square of the norm of the ground state 
wave function (denominator of the correlation function) was found in H, 



N N 

(o\U c (h)I[ b (h)\o) 

3=1 3=1 



JV 



>N>j>k>l 



N N 



II II HVjiVk) detiv 



U'=l k=l 



9( fj 

dfJ-k 



(2.21) 



Here dipj/dfi^ is N x N matrix 

dtpj_ 
dfJ-k 



Sjk 



N 



i=i 



(2.22) 



Let us emphasize that det(dipj/d/j,k) > (see Section 1.2 of Q). The thermodynamic limit of the 
square of the norm can be described by the following formula: 



lim 



det^ \ 



nf = i27rLp(^; 



del ( / —K 



(2.23) 



where K is an integral operator acting on some trial function /(A) as 



(Kf)(X) = J K(\,rif(n)dfi. 



(2.24) 



The proof can be found in S (see also Section X.4 of ||). 
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In order to calculate the correlation function we shall also need a description of excited states. 
We need to consider excited states which have one more particle than in the ground state 

N+l N+l 

J] B(Xj)\Q), and <0| J[ C(Xj), (2.25) 



3=1 

where Xj have to satisfy Bethe Equations 
a ( A i) j^ 1 /(Aj, Afc) = 1 



or 



5=1 



a (Aj) Tj 1 ^(Aj> ^fe) 



h(X k ,Xj) 



(-1) 



TV 



(2.26) 



d(A,-) /(Afc, Xj) 

k^j 

We shall further assume that the number of particles in the ground state N is even. In order to 



write the logarithmic form of the Bethe Equations it is convenient to introduce 



N+l 



(f>j = LXj + i In 



k = l 
k^j 



Xj — Afc + ic 
Xj - Afc - ic 



The Bethe equations can now be written as 



(2.27) 



(2.28) 



(fj = 2Trrij, 

where rij is an ordered set of different integer numbers rij+i > rij. One can prove that all Xj are 
real. In order to enumerate all the eigenstates in the sector with N + l particles we have to consider 
all sets of ordered integers rij . The square of the norm of the excited state is 

N+l N+l 

(oi n ^ n BiXjm = < 

3=1 3=1 



„N+1 



II 9{Xj,Xk)g{Xk,Xj) 

K N+l>j>k>l 



'N+l N+l 



n n K^j,h) det 



iV+1 



7=1 fc=l 



<9A fc 



(2.29) 



For the excited state det(d(pj/dXk) is also positive. We shall also mention that the scattering 
matrix of elementary excitations can be found in Section 1.4 of ||. It depends strongly on momenta, 
this shows that the model is not free fermionic. 



Now we can define the form factor in the finite volume 

N N+l 

f n = (oinco^M n B{Xj) 

3=1 3=1 

We shall calculate it in the next section. We shall also need the conjugated form factor 

N+l N 

<0| n C(A^t( M )]J 5(^)10) 

3=1 3=1 



(2.30) 



(2.31) 



e~ iht ■ exp 



JV 



(N+l 

*M£a|-£/4 



A' 



IX 



k=l 



I N+l _^ 
( E A 3 " E 



fc=l 
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Here we used the fact that the energy and momentum of the eigenstate are given by the expressions 

(2.32) 



N+l 

E N+1 = J2(^~h), 

3=1 



N+l 

Pn+i = Yj X r 

3=1 



(2.33) 



3 Form Factor 



The main purpose of the paper is to evaluate the correlation function. In the finite volume we shall 
use the notation 



t = <oinf =1 c^(o,o)^ ( x,t)nf =1 ^)io> 
(oinf=iC(^)nf=iB(^ 



(3-1) 



We shall use the standard representation of correlation function in terms of the form factors 



N N+l N+l N 

<oinc(MiMo,o) n b(\,)\o)(o\ n ^A^^n^oi ) 

E 3=l 3=1 3=1 3=1 
N+l N+l N N 

all {AW (0 | T] C(A .) T] B(A,.)|0><0| I] C(^) I] B (hM 



(3.2) 



3=1 3=1 3=1 3=1 

In order to calculate the form factor we need to know the action of the local field on the eigen- 
vector. This can be found in ||] (see also Section XII. 2 of ||). 



(3.3) 



N+l 


N+l 


1 ' N+l 


\ N+l 


mo) n £( A i)i°> = 


-i\/c Y a(Xe) 


n /(^,A m ) 




3=1 


i=i 


m— 1 


J m=l 



This permits us to represent form factor as follows 



N+l 



F N = -iy/c Y a (h 



N+l 



N+l 



m — 1 



m — 1 



N 



N+l 



x(oin^y n B(x m )\o). 

j=X m—1 



(3.4) 



Let us notice that the form factor is symmetric function of all the Xj because 

[B(X J ),B(X k )} = 0. 
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We now need to calculate the scalar product between the eigenvector and non-eigenvector 

N N+l 

(0| n C( Mj ) J] B(X m )\0), (3.5) 

j=l m=l 

where /i,- satisfy the Bethe equations, but X m do not. It can be done by the following theorem. 
Theorem 3.1 The following determinant representation holds for such scalar products: 

N N ( N 

<oi n c ^j) n = n ^m^) 

i=i i=i b=i 
II f(Aj, X k )g(fi k ,fj,j) >< Yl h(/J,j,X k )\ det(Mj k ) , (3.6) 

^N>j>k>l J [i,fe=l 



X 



where 



M ^ = g(ti k ,Xj) _ a(Aj) gjXjjjik) -p- /(Aj,/x m ) 
Jfc h(n k ,\j) d(Xj)h(Xj,fi k ) ^ 1 f(fi m ,X j y 



Here the spectral parameters {/J>j} satisfy the Bethe Ansatz equations (2. IS). The spectral parameters 
{Xj} are free and do not satisfy any equations. 

This theorem was proved in |1C]. 

For the scalar product, which appears in the expression for the form factor we get 

N N+l 

<oi n c(^) n BiXmm 

j' = l m=l 

N N+l 

= II ai^j^k)- n 9{Xk,Xj) -WW K^jiXm) 

N>j>k>l N+l>j>k>l j = l m=l (3.8) 

N N+l 
xj]^)' II d(X m )-det N M^. 

j' = l m = l 

Here the entries of the N x N matrix are 

M® H (\\*(\ ^ A /(Aj./Xm) j = l,... ,^-1,^+1,. ..,iV+l 

m=1 JlMm.) A jJ fc = l,...,iV 



Let us recall that 



x (*c) 2 , ,., a(A) 

' (A -» )= (A-,,)A-„ + < c) " nd = 



9 



Remember that the Bethe equations give: 



JV+1 



N 



™=n$£$ - n 



f(Xj, \X m , 
\ f(l*m, Xj] 



N 



(-i) N n 



h(Xj, [imj 



1 hiUm^j)' 



or equivalently, 



JV + 1 JV + 1 

°(^) II h (*ii^m) =d(\e) JJ h(\ m ,\ t ). 

m=l m=l 



Expression (|3.9| ) becomes 



=t(n k ,x j )-t(\ j , f i k ) I n 



JV+1 



/i(A p , A Jy 



JV 

n 



1 h{Xj,Xp) I \ rn=1 h(fj, m ,Xj) 



Using the obvious equality 

JV+1 

n 9(Xi,x m ) 

m=l 



II 9{Xe,X r 



\m=l 



JV+1 

II 9(Xe,X r 
K m=e+i 



/ l-l \ / JV+1 

II 5(A m ,A,) J] g(X e ,X 

\m=l / \m=e+l 



and substituting (3.8) into (3.4), we have 



JV+1 

FN = _ iV - c j2(-iy- 1 
t=i 



n 9(1*3, n 9(Xk,Xj] 

N>j>k>l N+l>j>k>l 



N+l N N+l 

x n h(x m ,x e )Yi n HvjiK 



m=l 

N 



7=1 m=l 



xjl^) II <A m )-det iV MW. 

j'=l m=l 



One can rewrite the determinant detAfM^ as 



/ 



det N M® 



N N+l 

n n 



\ j=1 h(n m ,Xj) 



N+l N+l 



p=l j=l 



where 



JV JV 

n h(nm,Xj) n h(Xj,fi r , 

s fk = Kl*k, A i)f^ *( A i' ^)^T 



II h(X p ,Xj) 

p=l 



n ^(A i; A p ) 

p=l 



j = l,...,£-M + l,. 

fc = l,...,JV 



10 



Let us substitute ( |3.15| ) into ( |3.14|) 



N+l N+l 

f n = -iy/c yi ai^j^k) n 9{hAj) n n 

N>j>k>l N+l>j>k>l m=l j=l 



X 



/AM-l \ AT JV+1 

]T (-l) /+1 det^ • n d (H) II d ( A -)- ( 3 - 17 ) 
\e=i / j=l m=l 



In order to simplify this expression let us study 

N+l 

Mi = Mi{\}=Y / (-l^detjvS^. (3.18) 
1=1 

Notice that Mi is an antisymmetric function of all {Xj} because Fn is symmetric and the product 
of functions g(Xk,Xj) is antisymmetric. In particular, 

Mi{X} = if Xj = A fc . (3.19) 

can be obtained from detS'^ +1 - ) by replacing A^ and Ajv+i- This is a special case of a 
permutation 

(Ai, • • • , A^, • • • , Xn, Aat + i) — > (Ai, • • • , Aat+i, • • • , Xn, Xg). (3.20) 
Since (— 1)^ _1 is the parity of this permutation, 

N 

Mi{X} = (-^IISPM 

Permutation of all {Ajy+i} J = l 

1 + det N (S jk - aS N +i,k)\ a=0 . (3.21) 



Here Sjk means from ( |3.16| ) and detSjk is the term & = N + 1 in ( 3.1? 
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— q- detN{Sjk - aS N +i,k) 



a=0 



is the sum of N terms where each of them differs from detSjk by the replacement of the l-th. line 
(corresponding to Xg) by the (N + l)-th line. We can use the expression ( |3.21 ) to simplify the form 
factor (|37T7|) 

N+l N+l 

f n = -%y/c Yi aiPjiVk) n 9^k,Xj) n n H^m,Xj) 

N>j>k>l N+l>j>k>l m=l j=l 

N N+l 

x JJ d(fij) JJ d(X m ) ■ Mi {A} . (3.22) 

j=l m=l 
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The complex conjugate of form factor is 

N+l 



N 



Fn = (0\ 1] 0(^(0,0) l[B( H )\0). 
Remember that c and all A, [i are real. Therefore complex conjugation gives 



ff(A,/x) = g(n, A), /(A,//) = /(/x, A), h(X,n) = h(/j,,X), 
t{X^ = t(fi,\), a(A) = d(A) = a _1 (A). 



So we have 



and for even TV 



>jfc 



Mi {A} = (-lJ^Mt {A} = Mi{\} . 
Hence for the complex conjugated form factor we get 

N+l N+l 

f n = iy/c yi g(fi k ,vj) n s{^j,^k) n n k^^™) 

N>j>k>l N+l>j>k>l m=l j=l 

AT AT+1 

x II ^) II a{\ rn )Mi{\}. 

j=l m=l 



For the correlation function the quantity \Fn\ is important: 

- ( N \ ( N+1 \ 

FnF n = c I JJ a(fj,j)d(fj,j) J y J| a(\ m )d(\ m )j 



( N 



2 = l,fc=l 



/ 



'AT+1 AT+1 



n 5(A,,A fc ) linn m^>^) i (-^m) 2 . 
i=i fc=i 



2 = l,fc=l 



or, 



FnFn 



N N N+l N+l 

<oi n c( Mj ) n B(Mi)io) • <oi n n bca^io) 
j=i j=i j=i j=i 



N+l N+l 



= C 



-2N \j =1 k=1 



n n h(Xj,X k ) ■ (Mi {\}Y 



INN \ g a - 

.n Jl HVj,Vk)J detAf^detTv+igg 
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This formula gives us |F/v| 2 &tx = t = 0. Also it is easy to "switch" on space and time dependence 
using the formula ( |2.31| ). Note that a(X)d(X) = 1 for nonlinear Schrodinger equation. Therefore the 
correlation function becomes 

Mo,o)v»t(M)>* = I -^r^—^— — E - Am l 3 ■ (3-29) 

(n n h^j,^)) -det N ^ {A}detiv + i^ 



3=1 k=l 



Here we used the new notation 



(N+l N+l \ (s^n+i , . ^JV , ,\ 

I] II KXsM {Mi{X}fe\^ (3.30) 
k=l 3=1 J 

where 

r(A) = itA 2 - iscA. (3.31) 



4 The idea of summation with respect to A 



Now let us consider the sum with respect to all {A}jy +1 in (|3.29|) 

^ A({A» 

w N+1 det N+i (§£) 

The idea of summation is the same as that which we used for impenetrable bosons (free fermionic 
case) ]ll| (see also Section XIII. 5 of ||). The factor (Aii {A}) 2 entering the r.h.s. of ( 3.30| ) contains 



(A + 1)! terms, all of them give the same contribution to the sum. So we can replace one of 
determinants Aii {A} by the product YljLi & 



03 

A({A}) 



E 

{x}n+i det N+1 



w+w e n e ~ rM ■ ( det ^lF 

f N+l N+l \ N+l / N \ 

n n MAi,A fe ) n eT(Aj) n ^ (4.i) 

k=l j=l J j=l \j=l J 

The sum with respect to all {A} Ar+1 means the sum with respect to all ordered sets of integers {rij} 
from Q2.28 ). We also can admit rij = because it leads to Xj = which does not contribute to 



j]) because of the antisymmetry of M.i{X\. The factor (A + 1)! is absorbed as 

N+l oo 

(A + 1)!]T = J] E • ( 4 - 2 ) 

{rij} i=l m=—oo 
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The correlation function becomes 

e -iht c -2N J| Ar _ 1 g-T&im) 



(^(0,0)^(x,t))jv 



(nf=infe=iMMi,Mfc))-det 



A' 



with 



'N+1N+1 \ N+l 

= inn h(\j,x k ) n ^ 

,i=i fe=i / j=i 

9 \ 

1 + — J detN(SjjSjk — aSjjSN+ik] 



(4.4) 



a=0 



The main difference between the free fermionic case (coupling constant c — > +oo) and the non- 
free fermionic case is that in the former it is possible to solve the Bethe equations ( |2.26| ) explicitly. 
On the contrary this is not possible for penetrable bosons. 

Our approach is based on the formula 

= £ / ^ +1 A A({A}) J] *(fc(A) - 2n nj ). 
Remember that detdtpj/dXk > 0. We shall also use the Poisson formula 

oo oo 

£ 6( x -2nn) = — £ e^. (4.5) 

n=—oo k=—oo 

So we have 

£ Lf ({A flU = E /^ +1 AA({A})n 1 ^(A)-2 7 rn J ) 

m---njv+i \°- e W+13A^/ j=l 



E (i) /^ +1 AA({A})n^ 



Thus we get the following representation for the correlation function 



-iht -2N yiN p -T(iim) 
dV™ -A\ llm=l c 



A- 



aii=infc=i%i,Mfe))-det^ 

/ n ^« n • («) 

ni— njv+i x 7 _ 3=1 \fc=l \ -" K > / 
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5 Quantum Dual Fields 



In this section we introduce the auxiliary Fock space and auxiliary Bose fields (f>o(X), <fii(X), 02(A) 
4>Aj (A) and 0n. (A) (j = 1,2). Further we shall call these operators dual fields |TI| (see also Section 
IX. 5 of ||). Dual fields help us to rewrite double products in terms of single products. 

By definition any operator (j) a (X) (a = 0, 1, 2, Ai, A2, D±, D2) is the sum of two operators: 
"momentum" p{X) and "coordinate" q(X). 

0o (A) = q (X) + Po (X); 

A .(A) = g^(A) +p Di (A); ^.(A) = g Dj (A) +p Aj (A); (5.1) 

<MA) = 91(A) + P2 (A); 2 (A) = g 2 (A) + Pl (A). 

All operators "momenta" p(X) annihilate the vacuum vector |0), all operators q(X) annihilate the 

dual vacuum (0| : 

P«(A)|0) = 0, (0| g a (A) = 0, for all a, (0|0) = 1. 
The only nonzero commutation relations are 

[po(A),<?o(a*)] = hi(h(\, fi)h(n, A)); 

[PD J (X),q D M] = S jk lnh(X, //); \p Aj (X),q Ak (ij)] = 5 jk ]nh(fJ,, A); (5.2) 

(We remind the reader that h(X,fi) = (A — /x + ic)/ic). The realization of these operators as linear 
combinations of canonical Bose fields is given in Appendix B. 

It is easy to check that all dual fields commute with each other 

[&,(A),&(m)] =0, 

where a, b run through the all possible indices. Using this property we can define functions of 
operators J-({e^ a ^}). One should understand such expression, for example as a power series over 
j e <MA) j_ 'p^g following simple formulas are useful: 



e PaW e q a (v) _ e ga(At)gPo(A) e [po(A),9o(A*)] ; (5 3) 

Mi Af 2 Mi M 2 

i=i fc=i j=l fc=i 
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M I M 



(Aj),g a (M)] 



10), 



A/ M / A/ 



(5.5) 
(5.6) 



Here {A}, {/i}, {/?}, {a} are arbitrary complex numbers, T is a function. One can easily prove these 
formulas. 

Let us define the very important dual field ip(X) as 



*/>(A) = ^o(A) + </> Al (A) + 0d 2 (A) + 02(A). 



(5.7) 



Theorem 5.1 The correlation function (4- r i) can be presented as the following vacuum expectation 
value in auxiliary Fock space 



(ii>(0,0)1>\x,t)) N 



iht„~2N N 



e c 



where 



and 



tlo1 V 7777f 



(0| Y[ (e po ^ m ^ e Pl ^ m ^ 



T^FTT E / ^ +1 A(7i(A^ +1 ) + 



9a 



x detAr (SjjSj k ^i(Xj)j2(lJ-j)l2(lJ-k) 
- aS jj S N+lk ^i(\j)ji(X N+1 )%(iJ,j)^2(^k)) _ n |0). 



a=0 



(5.8) 



(5.9) 



7i (A 3 -) = e iiA j™j+ T ( A j)+'/'( A j)+™^i( A j) 5 

Proof. Let us move factors 71(A) and 72 (aO ou t of the determinant in ( |5.8| ). In the r.h.s. of ( |5.8| ) 

we get 



1 + 



d_ 

da 



N 

(01 n [■ 



gP0Om) e Pl(Mm) 



m=l 

N+l N 
x J~J LikA m n m +T(A m )+V>(A m )+rim<MA m )j J~J j^ g — r(/Lt m )— 0(At m ) 

m=l m=l 



X 



detAr (SjjSjk - aSjj§N + ikj a _ |0) 



(5.10) 
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Using ( |5.4D we find 

N 



(oi n [ 



gP0(Mm)gPl(Mm) 



m=l 



N+l 



N 



m=l m=l 



|0) 



n n 



(5.11) 



Using fl5.5|) we obtain 

N+l N 

(Q| jj ^ e </>A 1 (A m ) + 0D 2 (Am)] 



detAr [SjjSjk — aSjjSN + ik) |0) 



m=l m=l 

= detAr {SjjSjk ~ aSjjS N+ i k ) ■ (5.12) 

Combining ( [5,11 ) and ( |5.12| ) we get the r.h.s. of ( |4.7| ) with A ({A}) defined in ( |4,4j ). The theorem is 
proved. 

Now we can rewrite the r.h.s. for (5JS) as follows 

oo 

* £ /^(ww + l) 



(2*r)' 



ni— njv+i. 



xdetiv (SjjSjkJiiXj^iVj^iVk) - aSj j S N+ i k %(Xj) A f 1 (X N+1 ) A f 2 (lJ'j) A f2{^k)) 

' oo 

i- 53 / dA Ar+1 7 1 (A iV+1 ) + ^ 



a=0 



«-iV+l_ 



(oo 
— 53 / dX j S j jSj k ii(Xj) A / 2 (lJ>j)72(tJ>k) 

i 00 

"°4^2 X! _/ ^A r +i^j^ii^A f +ifc7i( A i)7i(^A r +i)72(/ u i)72(Mfe) 



rij,njv+i_ 



(5.13) 



a=0 



In this formula we can perform the summation over integer {rij}. Indeed, rij enters only in 
function 71 (Xj): 

7! (A) = e iLX J n j+ T ( x i)+^( X j)+ n jM^j) ^ 
Recall that all dual fields commute with each other: 

[0a(A),&(/O]=O. 



This means that we can treat operators <j> a (X) as diagonal operators. Due to the formula (|B.8| ) from 
Appendix B we can consider operator i<^>i(A) as real function of A. Hence we can use formula ( [4.5D 
to sum up with respect to rij 



-1 00 00 

_ e in(L\-iM x )) = S(LX-ufn(X)-2irn). 



(5.14) 
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It means that A = A ra , where X n is a root of the equation 



LX n - 2im = i(/)i(\ n ). 



(5.15) 



The expression ( [5.15 ) is an operator equality, which is defined only on vectors of the form 
J] e <MA™)| ). Therefore one should understand this equation in the sense of mean value: 



(0\(LX n -iMK) - 2vrn)ne </ ' 2(Am) |0) = 0, 



(5.16) 



where {A m } are arbitrary real parameters. Then we can rewrite equation ( 5.16 ): 

A m — A n + ic 



An — X m + ic 



(5.17) 



h(X n , X m 

The r.h.s. of ( 5.17] ) is a real bounded function of A n . Moreover it is a decreasing function of A n , 
because 



d . v-^ . ( X m - X n + ic 



E 



2c 



< 0. 



5 An " ^ ~* V A n - X m + ic) ^ (A„ - A m ) 2 + c 2 
The l.h.s. of ( 5.17 ) is a linear increasing function of X n , hence equation ( |5.17|) has one real solution 
and this solution is unique. Also we have 

2c 



(0\(L - i^Xn)) n «^<*»>|0) = L + ^ tt _" 2 , > 0. 

m m \ n m) ~r~ 



Therefore, one can write 



*(LA-»fr(A)-2irn) = 



(5.18) 



(5.19) 



where A n is solution of equation ( |5,15 ), 
Later we shall use notation 

27rp(A) = l-^;(A). 
We now arrive at the following formula for the correlation function in finite volume 



(5.20) 



(V>(0,0)V> f (x,t)) N 



iht„-2N N 



e C 



(oi n 



m=l 







x ( G N (x,t) + — ) • det N (U jk - aQjQ k )\0) 



a=0 



where 



1 oo 1 
G N (x,t) = - ^7 



V(A„)+r(A„) 



L n^oo 2*-p(An) 



(5.21) 



(5.22) 



and 
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U» = 7 E 



V(An)+r(A„) -i(V'(Mj)+'/ , (Mfc)+T(Mj)+^(^fc)) 



L „f^oo 27rp(A„) 
x {t(Mfc,A n )e-^i^-t(A„,Mfc)e-^^} 

x {t( j u i! A n )e-^ 1 ^ -t(A n ,/ij)e-^(^)} , (5.23) 



1 00 p V>(An)+r(A„) -i(t/.(^)+r( Mj )) 



^ „±^oo 2vrp(A n ) 
x {t{^,X n )e^^ -t{\ n , N )e-^^)y (5 24) 



Formula ( ^.21 ) is the determinant representation for the quantum correlation function in a finite 
volume. 

6 Thermodynamic limit 

In order to calculate the correlation function in the ground state one should consider the limit 
where the number of particles and the length of the box tend to infinity with fixed constant density: 
N —* oo, L —* oo, N/L = D = const. In this limit the parameters {A n } are described by 
distribution density /3(A) 

^ -hi}- i 0;(A) ) • 

(see Appendix A). The sums in the expressions for Ujk and Qj can be replaced by corresponding 
integrals. Let us introduce the new function Z(A,/x) 

e -<t> Dl W 6 -<Pa 2 W 
Z X,fJ/ h(u,X) + h(X,u) ' 



Then we can rewrite ( |5,23D and ( |5.24 ) as 



1 ~ ( ic )2 e ^(X n )+T(X n ) 



L n f^ QO 2vrp(A n )(A n - Hj)(X n - Mfc) 
xe 4W(ft)+^f*)+ T fe)+-fe))z(A n)W )Z(A n , ft ), (6.2) 



1 ~ i Ce V( A nl + T(A n ) e ~^ 



V(A„)+T(A n ) p -i(V(M)+r(^)) 



Here Qj = —Q(pj), Qk = —Q{uk)- Using formula ( |A.25| ), we get 



icy2irp{pj) 
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(IC 



L 







2 ^ oot 2 f( ^ ) % 



(ic) 2 l-5 jk 
2 Hj - /ife 



; |(^(w)-^(^0+-fe)-r(M fe )) Z(M . )Mjfc)Z(/ ,. )/i .) cot ^| (/1 .) 



e 2 



\{'<l>{V-i)-'<l>{P-k)+T(jl3)-T(jl k )) 



Z(fJL k , fJLj)Z(ji k ,fJ, k ) cot —i(nk) 



+ 



(ic) 



2ir(fij - fj, k ) 



dX 



= ^(A)+r(A) 



A - /Uj A - /x fe 



xe -^W^)+T(^)+^ fc )+r(M fc )) Z ( Ai/Ufc ) Z ( Ai/ij ) + (1/L). 

Here we denote the principal value by the symbol 



A — n 



V.P. 



dX(-) 



1 

+ -x 



dA(-) 



rfA(-) _ 1 

\- H~ 2 J \- v + iO ' 2 J X — fji — iO 



Using ( |OlD we calculate the sum ( |6,3| ): 



Q(m) 



./" c?A 



27T / A — \i 



3 V(A)+r(A) e -i(V(M)+r( At )) z ^ Ai x 



~e^^+ T ^Z(fi, /i) cot + 0(1/L). 



Function £(//) is defined in Appendix A as 



(see (|A.2[ )), and hence 



L *■ 
cot -£(//) 



sin2 Triin) 



£0) =/*- yMv), 



e ii/i+^l(A«) _)_ J 
' e iLfjt+<t>i (fi) _ l ' 



2-e 



2^ vrv 4 

Let us turn back to the formula (|5.21 ), We can move all e Pl ^ m ^ to the right vacuum |0) 
each operator </>i entering into U and Q should be replaced by the rule (see ( |5.5| )) 



N N , / n AT 



m=l m 



n> 

m=l 



Taking into account Bethe equations (|2.1 



N 

n 



h(fJ, m ,IJ,j) 
j h(jJLj,Hm) 



T, (JV — even), 
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we get 



N 

11 e Pl ^ m) cot 

m=l 


= 


.W_(/Xj) 


N 
m=l 


= 


V 2 



n 



2 iV 



m=l 



pPl(Mm) 



where 



w± ( /U ) =e ^M/2 ±e -0lM/2_ 

Operator 27rp(A) also contains 0i(A), so it does not commute with pi(n): 



N 



N 



,Pl(Mm) 



where 



m=l m=l 



\ m=l / 

Hence we have the new representation for the correlation function 



e -iht c -2N 



detjv 



N , r, 

x(0| J] eP»(^) ( GN ( X}t ) + — 

m=l ^ 

xdet N (u jk - aQ(nj)Q(n k ) + 0(1/L 



a=0 



where 



Ujk — LSjk 2 



i ( ic ) 2 § "-(ft) 9 
2 ^ uj + (p,j) dfij 



Mfc —Mi 



—2 



,( ac) 2 1 - S jk 
2 /ij - y,k 



e 4»fe)-*(w)^fe)-(w))z( ft ,^.)2( ft)W ) 



w+(/ifc). 



27r(/i j - n k ) 



dX 



>(A)+r(A) 



A — Hj A — /ifc 



xe 4^(w)+ T (ft)+*(w)+^ft))z(A, Mfe)^(A, nj). 
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oo 

ic J dX 



2ir fX — ji 



= V(A)+r(A) e -i(VM+r(M)) 2 M \ 



2 v+W 



(6.10) 



Let us simplify the formula (|6.9|) and (|6.10|) . First, in (|6.9|) the term proportional to 1 — is 
defined only for j ^ k. Let us continue this term for all j and k using the l'Hopital's rule for j = k. 
Then 

{ic) 2 2np L (n j ) ry2l 



Ujk — L5jk~ i , 



.(ic) 2 1 
2 Hj - (i k 



I (^j )-*G* )+r(/* )) Z(/ij . , ^ ) Z(Atj . s ^. ) 



(ic) 2 



27T(/X j - /i fe ) 



= </>(A)+t(A) 



A — /ij A — /ijt 



where 



Using the Sokhodsky formula 



xe -^(ft)+^ 3 )+^te)+T(ft)) Z ( A) /i fc )Z(A, Mi), 



(6.11) 
(6.12) 



V.P.i = 7 — 777 i Z7T(^(x), 



x a; ± iO 

one can rewrite the expressions (lOll) and (Islob as follows 
(icf2irp L (nj 



Z (^,Mi) 



(ic) 5 



dX 



C.2 



|MA) 



+ 



e 2 



501(A) 



C.2 



5^1 (A) 



e 2 



501(A) 



27r(/Xj — /ifc) 7 ^+(A) V A — fij + iO X — \ij — iO X — p k + iO X — p k — iO 

— oo \ • • / 

Xe ^(A)+r(A) e -i(^ J O^(M,0-Hfr(Mfc)-K(Mfc))^(A j/ifc )Z(A,/i i ), 



(6.13) 



<ZA 



e ^i(A) e -|0i(A) 



+ 



; + (A) \A — ^ + iO A — /i — iO 



^(A)+r(A) e -l WW +rM) 2(A)/j)i (fu4) 



Now let us move the term proportional to the length of the box L out of determinant: 



N 



det N (U jk - aQ(fij)Q(nk)) = H 



a=l 



(ic) 2 2-Kp L (p a )L ( 



xdetA? 



V w+(Mo) 

- aQ(pj)Q(pk) _ a; + (//j)u; + (//fc) \ 
(ic) 2 2irp(p k )L Z(fj, j ,p, j )Z(fj, k ,fj, k ) J ' 



(6.15) 
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One can move the product Il^Li (ijff/ry) to the left vacuum (0| : 



N 



(°l IP'" 



a=l 



fa) ( Z (Va^a) 



(oi n 



X ( e -PA x (l*a) _|_ e -PD 2 (Va) \ 2 

a) 1 



oPo(P-a 



o=l 
N 



I e V2{Pa)/2 _|_ e -p 2 (Pa)/2 



(oi n^(^)- 



(6.16) 



0=1 



Now let us move the product in the r.h.s. of ( |6.16| ) to the right vacuum. The only operator in the 
expressions for U and Q which does not commute with V{p a ) is = </>o(A) + (A) + 0d 2 (A) + 

N 

02(A). In order to move Y\ V(fi a ) through the determinant we use the following lemma. 

o=l 

Lemma 6.1 For arbitrary M = 1, 2, . . . and arbitrary complex numbers Ai . . . , Xm, (3\ . . . , /3mV 



m 



M 



m=l 



W(X m ) 



|0) 



m=l 



Proof. The proof is straightforward: 



Villa) n 



m=l 



I e P2(p.a)/2 _|_ e ~P2(p.a 



)/2 



n 

m=l 



|0) 



(6.17) 



M M 

JJ e^( A -)|0) J] [/i(/"a,A m )/ l (A m ,/z a )]^ 

m=l m=l 



/ M M \ 

' II WWa)]""*^ II [/l(/"a,A m )]-^ \ 
m=l m=l 



Ai 
\m=l 



+ n 

m=l 



h(flg,Xm) 
h(\ m ,Ha) 



P m /2 



J 



M 
m=l 



|0). 



This proves the lemma. 

Since the determinant in r.h.s. of ( 6.15| ), being a function of operator ip, is some linear combina- 

M 

tion of products of the type n e 13 " 1 ^" 1 ^ (with different M, {/3} and {A}), we can move ]la=i V(p a ) 

m=l 

to the right vacuum without changing the matrix elements of the determinant ( |6.15| ). Therefore we 
have 

„—iht 



N 



U^npML) (G N (x,t) + —) -det N (W jk + 0(l/L 2 )) |0) 



0=1 



a=0 



where 



1 



(V jk - aP{nj)P{n k )), 



(6.18) 



(6.19) 
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and 



V, 



jk 



dX ( esfcW e-^W 



w+(A) 



+ 



E2 



e 2 



30i(A) 



X — fij + iO X — Hj — iO X — Hk + iO X — pLj^ — iO J 



(6.20) 



^0) 



e 2 



3</>i(A) 



27rZ(/x, /i) 7 ^+(A) I A — ;U + iO ' X — fj, — iO ^ 

Xe ^(A)+r(A) e -i(^(M)+r( M )) z(A)/i) _ 



Recall that in the thermodynamic limit 



det 



A' 



iV 



n (27rp0i o )L) det(J - —K), (see (2.23)) 



0=1 



G N {x,t) -> G(x,t) = ^ / e^M^M di/, 



(6.21) 



(6.22) 



(6.23) 



px(A)-p(A), (see©), 

The formula ( |6.18| ) contains the determinant of the matrix W . In the thermodynamic limit it 
will turn into a determinant of an integral operator. The simplest way to see this is to express det W 
in terms of traces of powers of the matrix {V — aPP). The replacement of summation by integration 
(in the limit) is straightforward and is explained in detail in Section XI. 4 of |J. Therefore the 
determinant tends to the Fredholm determinant. Now we arrive at the main theorem. 

Theorem 6.2 In the thermodynamic limit, the time- dependent correlation function have the 
following Fredholm determinant formula. 



/ d \ det (i + ±-V a 

(V(0,0)^(^t)) = e-^(0| (G(x,t) + — X 

V da J act (l-^-lT 

Here the integral operator V a is given by 



(6.24) 



a=0 



(V a f) (A) = f (y(A, M ) - aP(ti)P(\j) f{fi)d(i, 



(6.25) 
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where q is the value of spectral parameter on the Fermi surface. Here the kernels V (111,112) and P(fi) 
are given by 



and 



where 



00 



x 



2n((ii - fjL 2 )Z(m,iJLi)Z(n2,H2) 
fl\ ( e i<M A ) e -|^W e ^i(A) e -^i(A) 



w+(A) \A — /xi+iO A — /ii — iO X — fi2 + iO X — [12 — iO I 

x ^(A)+r(A) e -I(V>(Mi)+r(Mi)+^(^)+r(M2)) Z ( A) ^) Z (X, fix), (6.26) 



A , . 7 dX f ei^W e-3^W 



2irZ([x,n) J oj + (X) yX — fi + iO X — fi — iO I 

Xe *W+rW e 4^M+^))2(A, M ) ) (6.27) 



w+ ( j u) = e^W/ 2 +e-^W/ 2 , 



/i(/u,A) /i(A,/i) 
r(A) = iiA 2 — zccA. 

The integral operator K is given in (2.24). 



We want to emphasize that formula ( |6.24 ) is our main result. 



It is easy to show that it has the correct free fermionic limit. He—* +00 (free fermionic case) 



then all commutators (5.2) of auxiliary "momenta" and "coordinates" go to zero because in this 



limit h(X, fi) — ► 1. Hence one can put all dual fields </> a (A) = 0. In particular 



V>(A) = 0, 0i(A) = O, 
w+(A) = 2, Z(X,fi) = 2. 

whereby we have 



V(»i,»2) C =°° - ( ' , fd\(-^—- \~~ — ) e«X>-ftrt-fr>»\ (6.28) 
7t(a*i - H2) J \ A - m A - ii 2 J 



-DC 



00 

P( M )^-I / dA ^_ e -(A)-|rM. (6.29) 

7T J X — 11 
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K = 0, (6.30) 

oo 

G(x,t) c =°°^- / due T{u) . (6.31) 

Z7T J 



Substitution of these formulae into ( |6.24 ) reproduces the result of 11] - In order to obtain Lenard's 



determinant formula M one should also put t = 0. 



Summary 



The main result of the paper is formula (6.24). It represents the correlation function of local fields 
(in the infinite volume) as a mean value of a determinant of a Predholm integral operator. In order 
to obtain this formula we introduced an auxiliary Fock space and auxiliary Bose fields (all of them 
belong to the same Abelian sub-algebra). This is the first step in description of the correlation 
function. In forthcoming publications we shall describe the Fredholm determinant by a completely 
integrable integro-differential equation. Then we shall solve this equation by means of the Riemann- 
Hilbert problem and evaluate its long-distance asymptotic. 
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A Summation of singular expressions 

Let us consider equation ( |5,15|) 

L\ n - iMK) = 2nn. (A.l) 
where i</>i(A) is a real and bounded function for ImA = 0. Let us introduce a function £(A) 

£(A) = A-i0i(A). (A.2) 

Obviously 

, . 2irn 

£(An) = — • (A.3) 
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Comparing with (5.20) we get 
It follows from equation (A.l) that 



where 



2vrp(A) = l- I 0' 1 (A)=e / (A). 



|A n+1 -A n | < -(tt + M), 



M= sup \(f>i(\)\- 

— oo<A<oo 



(A.4) 



Hence, |A n +i — A n | — > if L — > oo and we can make the following estimate 

' /5(A n ) + 0(l/L 2 ). 



L{\ n +i — A n ) 



(A.5) 



Due to ( |5.18[ ) we have 27r/5(A) > 0. 

During study of thermodynamic limit the following sums appeared 



1 



S =T E 



L n f^ oo 2vrp(A n )(A n -//)' 



(A.6) 



Here /(A) is some smooth function, \i is some fixed point on the real axis. We shall be interested in 
the asymptotic of this sum when L goes to infinity. 
Let us present ( |A.6| ) as the sum of three summands 

Ni-1 



1 



E 



/(•V 



2vrL \ fr^ /3(An)(A„ - /u) 



+ 



/(Ar 



+ y ^ + y 



(A.7) 



Here iVi and are integers such that in the limit L — > oo, the following properties are valid 

< /x — Ajv-l < oo, < Aj\r 2 — M < 00 • 



(A. 



Obviously the first and the third summands in (A.7) have no singularities in the domain of summa- 
tion. The corresponding sums are integral sums, for example 



Si 



L 



Ni-l 

E 



/(A») 



2vrp(A n )(A n - fi) 



j y ( A„)(w-A„) +0(1/L2) 

— oo 



(AS 
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An analogous formula is valid for S3 



E 



/(An) 



L n=A^+l 27r /5(An)(A n ~ /i) 2?T J A-// 



/(A) 



dA + 0(l/L). 



(A.10) 



iV 2 



Consider the second summand in (A. 7) 



So 



N 2 

E 



L „=jvi 27T/3(A„)(A n -//) 



One can present S2 in the following form 



C _ C-(l) , c(2) 
J2 — «->2 + '-'2 ) 



where 



5 (D = 1 g / /(An) /(At) 



£ n t^ V 27r P(An)(An " M) £(A n ) " 

s (2) = M v 1 (a id 

£ n^^An)-^)' 

Due to Q 5^ has no singularities in the domain of summation. Therefore it can be replaced by 
the corresponding integral 

s m = ±_ I ( m _ M»/M dA + (A 12) 

2 V A -M f(A)-«M/ 

Using ( |A.3| ) we can rewrite in the following form 

s (2) = ^ 1 



2?r n = Ni n-±^) 

The last sum can be calculated explicitly in terms of the logarithmic derivative of the T-function: 

ip(x) = — lnT(x). 
ax 

We shall use the following properties of the ^-function 

ij)(x) + - = il>(x + l), (A.13) 
x 

lp(x) — lp(l — x) = — IT COt 7TX, (A. 14) 

ip(x) -> lnx + 0(l/x), x — > +00. (A.15) 
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Now using ( |A.13 ) we can write 



S. 



(2) _ f(fj) 



2vr 



^ N2 _ L.fo) + 1) _ _ 



(A.16) 



The argument of the second ^-function in ( |A.16| ) is negative. Using ( A.14| ) one can flip the sign 
of this argument 



S. 



(2) _ f(fi) 



2tt 



^( N 2 - ^10-0 + 1) 



L 



V^£G«)-iVi + l)-7r cot 



'2vr 



L 



(A.17) 



Remember now that < fi — Xn ± < oo and < Xn 2 — fi < oo (see (A. 8)). This means that the 



arguments of ^-functions in (A.17) tend to infinity if L — > oo. Therefore we can use the asymptotic 
formula ( A.15| ) 



,(2) _ Hm) 



2tt 



In 



+ 0(1/L). 



(A.18) 



Now let us turn back to ( A. 12] ), Let us present the r.h.s. as the difference of two integrals. Both 
of them should be understood in the sense of principal value (V.P.): 



c(i) _ ' 

<3 2 — ~ 



1 



/(A) 

2-7T f X — fj, 



\N 2 



dX 



C(a)-I(m) 



dA + C7(l/L). 



Due to ( [OD one can compute the second term in ( |A.19| ) explicitly 



(A.19) 



dA 



/(M)/5(A) /(/,) 



£(A)-£( M ) 2vr /£(A)-£( M ) 



de(A) 



Combining now QA.18; ), ( |A.19[ ) and ( |A.20Q we get 



\]V 2 



Finally, using ( |A.9 ) and ( A.lOj ) we find 



S 



L . 27r / 5(A n )(A n - /u) 



E 

=— c 

2vr /. 



/(A, 



n=— oo 
oo 



A _ dA-Moot||(M) + 0(l/L). 



(A.20) 



(A.21) 
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This formula describes the asymptotic behavior of the sum ( |A.6| ). 

For the evaluation of the thermodynamic limit of our determinant representation it is necessary 
to consider a sum containing the second order pole 



-. oo 



/(An) 



L n fr' 00 2vr/5(A n )(A n - n) 2 ' 



(A.22) 



Taking the derivative of ( |A.21| ) with respect to fi, we get 



5' 



L 



4sin 2 lf£(/x) 



(A.23) 



We can use formulae ( A.21 ) and ( |A,23 ) to calculate thermodynamic limit of (|6 



-. oo 

- y 



f(X n \fij,fi k ) 



2np{\ n )(\ n - Hj)(\ n ~ Mfc) 



1 



+ 



2(Mj - Mfc) 
1 



f(Hj\Hj,fJ, k ) COt - f(fJ, k \fJ,j,Hk) COt ^ri{Hh] 



d\f{\\nj,n k ) 



1 



1 



+ C(1/L). 



(A.24) 



27r(/x.,- - /j, k ) f "' v " Irv ' r ~' v/ V A - fj,j A - [i k 

— oo x ' 

One should understand the r.h.s. of this equality by l'Hopital's rule if j = k. It is also useful to 
extract explicitly the term proportional to the length of the box L 



-. oo 

- y 



2vrp(A n )(A n - fij)(Xn - Hk) 



8j k L 



2irp(n j )f(fj, j \fi j ,n j ) 
4sin 2 -§^) 



27r(/Uj 

1 - Sjk 
2 (p-j - Mfc) 



f(jij\fJ,j,fMk) cot -^iifij) ~ f(Hk\Vj,Vk) cot ^i(nk) 



-f cot^)^- 
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f{}ij\lij,Hk) - f(Vk\Vj,Hk) 
We have used formula? ( [A. 21 ) and ( |A.25| ) in Section 6. 



+ 0(1/L). 



(A.25) 



B Representation of dual fields 

What is the relation between our dual fields and the canonical Bose fields. Canonical Bose fields 
ipi(X) can be characterized as follows 

[^(A),^(A)]=W(A- / u), (B.l) 
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(do not confuse this ^ty(A) with the dual field ip(\)) and 

^(A)|0) = 0, (0|V/(A) = 0. (B.2) 

The dual fields which appeared in this paper have the form 

0a(A) =q a (X)+Pa(X), (B.3) 

where p a (A) is the annihilation part of ^ a (A) and q a (X) is its creation part. Their commutation 
relations are 

[Pa(A), qb(fi)} = «afe(A, //). (B.4) 
Here a a fc(A, n) is some complex function. 

p a (A)|0)=0, (0| g a (A) = (B.5) 

One can represent our p a and qt> in terms of ipi and ifij, for example as 

Pa (A) = Va(A), 

(B.6) 

oo 

#>(m) = J2 dva ch {y,[i)^\{y). (B.7) 

c 

— oo 

This shows that the dual fields, which appear in this paper are linear combinations of the standard 
Bose fields. 

Let us now consider a related issue. We can realize the dual fields (fii(\) and 4>2 (A) as 

q 2 (X) = 4(X), Pl (A) = Vi(A); 

oo oo 
— oo — oo 

Here f means Hermitian conjugation, and 

[Vi(A),40u)] = [^ 2 (A),^(/x)] = 5(A-/x). 

Other commutators are equal to zero. These commutation relations differ from flB.l| ) only by a 
trivial relabeling. Then 

oo 

<MA)=V>!(A)+Vi(A), MV= J hx^^(4(y) + th(y))du. 
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This means that 2 (A) and i0i(A) are Hermitian operators: 
4(A) = 02(A), (#x(A)) f =#i(A), 
After diagonalization they will turn into real functions. 



for Im A = 



(B. 



C Reduction of number of dual fields 

We would like to reduce the number of dual fields in the determinant formula for the correlation 
function in ( |6.24 ). Here we shall show that 



and 



2 (A) = 0. 



Recall the definition of the dual quantum fields fl5,l|) and (5.2) 



0o (A) = qo(\) + po{\); 

0A ] (X)=q A] (X)+p Dj (X); 

0i(A) = ?1 (A)+p 2 (A); 



0^(A) = ^(A)+^.(A); 
02(A) = g 2 (A) +pi(A). 



(C.l) 



(C.2) 



[Po(A),go(/-0] = hi(h(\,n)h(n,X)); 
[PD 3 {X),q Dk {^)} = S jk ln h(X,fi); 



[Pi(A),gi(/i)] = In 



\PA 3 (X),qA k (v)} = 5 jk lnh(fj,,X); 

h(n,X) 



h(ji,xy 



[P2(A),g 2 (M)] = In 



h(X,n)' 



Remember also the definition of the field V'(A) (5.7) 

^(A) = O (A) + Al (A) + D2 (A) + 2 (A). 

Notice that q\(X) and p 2 (A) entering into 0i(A) do not commute only with ip(n): 

h(n,X) 



Mv),qiW] 



hi 



In 



/i(A,^)' 
h(fi,X) 



(C.3) 



h{X,n) 

On the other hand qD x {X) — <M 2 (A) and pa 1 (X) — pr> 2 (A) entering into 0_D!(A) — 0^ 2 (A) also do not 
commute only with V'(/ u ) : 

h{n,\) 



[( P A 1 (X)- PD2 (X)),i;(fi)) 
[^),(q Dl (X)-q A2 (X))} 



hi 



hi 



h(X,|S) , 
h(X,fj,) ' 
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so we can identify 

1 (A) = ^ 1 (A)-^ 2 (A). (C.4) 

Then we can also put ^2 (A) = 92(A) + pi(X) = because after the replacement (|C,4| ), operators 
92(A) and pi(A) commute with everything. 
Such a replacement implies 

o, + (A) = e5^i( A )+^2( A ))z(A,A). (C.5) 

It means, that Fredholm determinant in ( |6.24D realy depends on three dual fields i/j(X), (f>A 2 {X) 
and 0o 1 (A). We shall use this fact in our next publications. 



D Thermodynamics 

The thermodynamics of the quantum nonlinear Schrodinger equation was described by C. N. Yang 
and C. P. Yang p3| . It involves few equations. The central equation is for an energy of excitation 
e(A): 

00 

4X) _ A * _ h _ I j ^-l—, ln (, + .-*) ^ (d.d 

—00 

Other important functions are the local density (in momentum space) of particles p p (X) and the 
total local density pt(X) (it includes particles and holes). They satisfy equations: 

00 

2vr Pt (A) = l+ J — *L-— Pp fr)d fi , (D.2) 

—00 

^-^l + e^V^^A). (D.3) 



Pt(X) 

The global density D = N/L can be represented as 



00 



D = / p p (X) dX. (D.4) 



-00 



In order to obtain the determinant representation of the temperature correlation function we can 
use the following representation 

wo , w t (x , t))T = tr(«-»WQW'(».0) = WM)W, (D . 5) 



tr e t 



(n T \nr) 
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Here |0y) is one of eigenfunctions of the Hamiltonian, which is present in the state of thermo- 
equilibrium. It is proven in Section 1.8 of Q that the r.h.s. of ( p.5| ) does not depend on the 
particular choice of 

Now we have to recalculate thermodynamic limit. First we shall return to the determinant 
representation of the correlation function in a finite volume, see formulas (6.18)-( 6.21| ). We should 
also notice that the thermodynamic limit of square of the norm should be changed (comparing to 
as follows: 

detjv^ - J] (2iTpML) det(J - —Kt), (D.6) 



0=1 



see Section X.4 of ||. Here pj correspond to \£It)- I n the thermodynamic limit summation with 



respect to indicis j and k in ( 5.18 ) will be replaced by integration p p (ph) dpk- Also Pl(/x) defined in 
(|6.12| ) goes to p t : 



After dividing p p (pk) dpk by pt(pk) which appears in the denominator (6.19) we shall obtain an 



integration J #(A) d\(-) insted of / dA(-). The details of these calculations are explained in Section 

— 00 — q 

XI.5 of II. 



The integral operator Kt can be defined by its kernel 



(D.7) 



,c 2 + (pi - p 2 ) 2 , 

Now let us formulate the final formula for the representation of the temperature correlation 
function of local fields in the thermodynamic limit 



(m 0)tf Or, *)> = e~ iht (0\ G(x, t) + — x 



d\ det(/+^T> a>T 



da ' det (l-±K 
Here the integral operator V a ,T is given by 

(v Q>T f) (A) = f° (v T (\,p) - aP T (p)P T (X)) f(p)dp. 



10) 



(D. 



a=0 



Here the kernel of ^ I Vt(//i, p-z) — aPT{pi)pT{p2) ) differs from the zero temperature case by the 
measure and limits of integration: 



VHMl, M2) - aPTGui)P T M = f^(Mi,A*2) - aP(p 1 )P(p 2 )) ^{^)^fd{p^). 



(D.10) 



It acts on the whole real axis —00 < p < 00. Here V(p\, p 2 ) is given exactly by ( |6.26| ), P(p) is given 
by formula (6.27) and G(x,t) is given by (6.23). 
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